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1 Introduction 


At thirty meters above the surface of Mars, a one-kilogram, twin-rotorcraft 
helicopter Martian scout (eyes in the sky) hovers without batteries. Its rotor 
engine is powered by a solar panel that is energized from an infrared laser 
mounted on the planet’s surface. For this simulation, the scout is restricted to 
vertical motion. Its mission is to use an onboard camera to seek out safe and 
interesting locations in which to direct a land rover (boots on the ground). 


Initially, with the solar panel illuminated, the scout hovers at an altitude of 
thirty meters. At some point in time, the decision is made to it return the craft 
to the surface. This is accomplished by momentarily blocking the laser beam 
until the scout drops to an altitude of about twenty-seven meters with a speed 
of four meters per second (see page 4 for calculations). 


When the scout reaches an altitude of twenty-seven meters, the beam is 
unblocked as the probe continues its descent with its solar panel fully 
energized. Note that a safe touch down speed is approximately one meter per 
second. The blocked-beam descent phase will be solved analytically to 
provide updated initial conditions. 


It is assumed that the laser power delivered to the solar panel is not 
significantly attenuated as it propagates through the Martian atmosphere. 
Thus, the solar panel output will not change and power to the rotors will 
produce a constant thrust until the scout touches down, at which time the 
laser would be deactivated. 


For this application note, quadratic drag is assumed. 


For comparison, analog computer, numerical, and analytical methods will be 
employed with the results displayed in table 1. 


2 Mathematical Modeling 


Phase 1: Blocked-beam descent (analytical method only) 


Starting with Newton’s second law of motion, and assuming down to be 
positive (forward), 


Ma = Fret 
Ma = —drag + weight 
Ma = -kv? + Mg 
a=-k/M v?+ g 
dv/dt =-k/M v?+¢ 
dv/dt = —-k/M (v2 - Mg/k) (1A) 
Setting dv/dt =0 and solving for terminal speed, vr, 
v4 = Mg/k (2A) 
Rewriting (1A), 
dv/dt = -k/M (v2 - v4) 
Using the chain rule for differentiation, dv/dt = v dv/dz. Thus, 
vdv/dz = -k/M (v2 - v4 ) 


where, z (assumed to be positive) is the distance traveled from its initial 
hovering position. 


Continuing, 


vdv/(v2 - v4) = -k/M dz 
V 


| vavsev? —va) = -k/M | dz (3A) 


0 
By inspection, the the right-hand side of (3A) becomes —kz/M. 

Solving the integral on the left-hand side of (3A) requires a little more effort! 
Setting p = v2 - v4, 

dp = 2vdv, or vdv = dp/2, and p(0) = 02 —v4 = —v4 


Substituting dp/2 for vdv, 


dp 
yz | 2 =1/2In (p/p 
p(0) 

= 1/2 In|(v? — vz) /(-vp)| 

= 1/2 In|1 - (v/vr)?| since v < vr, 

= 1/2 In(1 - (v/vr)?) 
Equating this to the right-hand side of (3A), 

In(1 — (v/vr)?) = —2kz/M 
Following a bit of algebra, 
v = vrV(1— e4(-2zk/M )) 

Letting k = 0.10 kg/m, M = 1.00 kg, g = 3.71 m/s?, and using (2A), 


vr = 6.09 m/s 


and 
v = 6.09 m/s V(1— e*(—0.20 m* z) 


At z = 2.82 m, v = 6.09 m/s V(1- e4(-0.20 m7 x 2.82 m)) = 4.00 m/s and 
altitude = 30 m - 2.82 m = 27.18 m = 27 m. 
Phase 2: Unblocked-beam descent 


Starting with Newton’s second law of motion, and assuming down to be 
positive (forward), 


Ma = Fret 
Ma = —drag + weight — thrust 
Ma = —-kv? + Mg-T 
a=-k/M v?+ (Mg-—T)/M 
dv/dt = —k/M v? + (Mg - T)/M 
dv/dt = —-k/M (v2 - (Mg —-T)/k) (1B) 
Setting dv/dt =0 and solving for the terminal speed, vr, 
vt = (Mg - T)/k (2B) 
Rewriting (1B), 
dv/dt = —-k/M (v2 - v2) 
Using the chain rule from calculus, 
vdv/dz = —-k/M (v2 - v4 ) 


vdv/(v2 - v4) = -k/M dz 


For mathematical convenience, z will be set from 2.82 m to 0 m. 


Vv 


| vav/e? —va)= -1yM | dz (3B) 
v(0) 0 


Following the mathematical techniques used in phase 1, 
the integral on the right becomes —kz/M and 
the integral on the left becomes 
1/2 In|(v2 — v2) /(v2() — v4)|. However, this time v(0) = 4.00 m/s.. 
Since v(0) and v > vr, 
1/2 In((v? — vt)/(v2(0) — vp) = —kz/M 
Following a bit of algebra, 
v = V((v?(0) — v4)e*(—2kz/M) + vz) 
Substituting numerical values, 
v= vV(15e*(—0.20 m!z) + 1) m/s 
Requiring vr = 1.00 m/s and using (2B), 
thrust, T = 1.00 kg x 3.71 m/s? - 0.100 kg/m x (1.00 m/s)?= 3.61 N. 
This is the combined thrust for the rotors. 
Substituting values for k, M, and vr, and omitting units for clarity, 
dv/dz = —0.100/1.00 (v2 - 12)/v 


dv/dz = —0.100(v2 - 1)/v 


3a Computer setup (patch cord version) 
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Figure 1: Computer setup 


























3b Computer setup (IC/discrete component version) 
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Red arrows connect to VCC. 
Blue arrows connect to VEE. 7411 





Green arrows connect to ground. 
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Figure 2: Basic breadboard layout for divider circuit 
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Figure 3: Basic breadboard layout for hovering on Mars 





4 Numerical Method 


(Euler’s Predictor-Corrector method using a hand-held programmable 
calculator) 


Code: TI-BASIC 


PROGRAM:HOVERING 
:-ClrHome:ClrDraw 
:’DV/DT=—0.1(V2-1)/V” 
“INITIAL CONDITION:” 
:(0-Z:4-V 

:”’STEP SIZE:” 

:0.5-H 

‘Fix 2 

:Lbl 1 

‘If Z>30:Then 

:Goto 2:Else 

:Disp {Z,V} 
:—0.1(V2-1)/VF 
‘Z+HZ 

:V+HF>W 
:-0.1(W2-1)/VS 
((F+S)/2A 
:V+HA-V 

:Pause 

:Gotol 

‘Lb! 2 

‘End 


5 Results 






































z (m) Analog computer Analytical Euler’s Predictor- 
from oscilloscope v (m/s) v (m/s) Corrector Method 
scale: 1s=1m estimated from v (m/s) 
oscilloscope Calculator 

00.0 4.0 4.00 4.00 

05.0 2.5 2.55 2.55 

10.0 1.8 1.74 1.74 

15.0 1.5 1.32 1.32 

20.0 1.3 113 1.13 

25.0 12 1.05 1.05 

27.2 touch down 12 1.03 1.03 





Table 1: Solution comparisons for unblocked-beam descent 
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Figure 4: Martian scout descent speed during powered phase* 


*For this application note, the oscilloscope display was produced during a 
single run using a differential equation analog computer constructed from 
operational amplifiers, analog multipliers, and discrete components with 

tolerances within 20%. 
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Figure 5: Basic breadboard layout 
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